Inverse cascade in Charney-Hasegawa-Mima turbulence 
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The inverse energy cascade in Charney-Hasegawa-Mima turbulence is investigated. Kolmogorov 
law for the third order velocity structure function is shown to be independent on the Rossby number, 
at variance with the energy spectrum, as shown by high resolution direct numerical simulations. In 
the asymptotic limit of strong rotation, coherent vortices are observed to form at a dynamical scale 
which slowly grows with time. These vortices form an almost quenched pattern and induce strong 
deviation form Gaussianity in the velocity field. 



The existence of an inverse cascade is the most remark- 
able property of two dimensional turbulence. It was pre- 
dicted by Kraichnan M for Navier-Stokes equation: as a 
consequence of inviscid enstrophy conservation, energy is 
forced to flow to large scales. The inverse cascade can be 
sustained only in presence of an external forcing inject- 
ing energy at a characteristic scale into the system. At 
scales larger than the forcing the turbulent flow is essen- 
tially random with Gaussian velocity difference statistics 
following Kolmogorov scaling 0. Thus, the presence of 
the inverse energy cascade prevents the formation of the 
large scale coherent structures observed in the case of 
decaying turbulence || . 

Large scale coherent structures in presence of the in- 
verse cascade have been observed only in presence of 
a characteristic scale breaking scale invariance. A well 
known example is the so-called Bose-Einstein condensa- 
tion, when the energy accumulates at the largest available 
scale j| forming vortices at the system size. 

Another example of vortex formation is the quasi- 
crystallization phenomenon observed in Charney- 
Hasegawa-Mima (CHM) turbulence, a paradigm for both 
geostrophic motion in planetary atmospheres (^] and 
drift-wave turbulence in magnetically confined plasma 
||. For the stream function ip(x,t) CHM equation is 
written as: 



0_ 

di 



(VV - A 2 V0 + J(V 2 ?/>, ip) = D + F 



(1) 



where J denotes the Jacobian and D and F are damp- 
ing and forcing respectively. In this case, vortices have 
been observed to form, in a quasi-crystal structure, at 
the intrinsic scale 1/A, corresponding to the Rossby de- 
formation radius in the atmosphere or to the effective 
ion Larmor radius in plasma H . 

In this Letter we focus on dynamics on scales much 
larger than A -1 . In this regime there is no intrinsic scale 
involved in the evolution of the system, which neverthe- 
less exhibits formation of strong vortices. In this case 
we observe that the scale of vortices is a dynamical one 
which increases in time as a consequence of vortex merg- 
ing. The characteristic time of evolution slows down lead- 



ing, in the limit of large Reynolds numbers, to a disor- 
dered pattern of quenched vortices. Despite the pres- 
ence of strong vortices, we find that the two-dimensional 
3/2 Kolmogorov law for the third-order velocity struc- 
ture function holds, independently on the value of A. As 
a consequence, the kinetic energy spectrum follows Kol- 
mogorov scaling but with a different constant with re- 
spect to the Navier-Stokes turbulence. 

The CHM equation ([!]) has two quadratic inviscid in- 
variants corresponding to total energy 



E = E k +E x = ^{(ViP) 2 + \ 2 ^ 2 ) 



(2) 



where (...) denotes spatial average, and total enstrophy 



Z = Z k + Z x = U(V 2 ^) 2 + \ 2 (V^) 2 ) 



(3) 



Both the inviscid invariants consist of two terms, the first 
corresponding to kinetic contribution and the second to 
potential one. The kinetic terms are, by definition, the 
only ones which survive in the Navier-Stokes limit A — > 0. 

The range of scales are separated by the characteris- 
tic wavenumber A. For k 3> A the kinetic contributions 
dominate in (|[^), at very large scales k <C A the leading 
terms are the potential ones. In the following we will 
assume that the forcing F is limited to a narrow band of 
wavenumber around kf. This will be the other relevant 
wavenumber in our problem. 

Kolmogorov-like dimensional analysis can be easily ex- 
tended to the present problem ]q-|lQ||. If kf 3> A we re- 
cover the well-known energy spectra for two-dimensional 
Navier-Stokes turbulence with energy spectrum E{k) oc 
k~ 5 / 3 and E{k) oc k~ 3 for A < k < kf and k > kf 
respectively. When kf <C A one obtains the predic- 
tion E(k) cx fc~ n / 3 for k < kf and E(k) oc fc~ 5 for 
kf > k > A. Dimcnsionally predicted spectra have been 
confirmed by direct numerical simulations p,0,0| but 
little is known about structure functions and probability 
density functions. 

The starting point for a statistical approach to turbu- 
lent cascade is the energy flux, written in the physical 
space as [O 



1 



W>(x + £) + A 2 V(x)V-(x + £))\ NL 

(4) 



where the subscript NL stands for the nonlinear contri- 
bution of (|l|) to the time derivative. Making use of (|l|) 
and of integration by parts we easily obtain 



e{£) = -<V(x) J(^(x + £), VV(x + £))) + (x 



(5) 



in which the parameter A has formally disappeared. As 
a consequence we can make use of the well known result 
for Navier-Stokes (i.e. A = 0) and write E3 



((Su L f) = -el 



(6) 



where represents the longitudinal increment of the 
velocity v = (d y ip, —d x ip) and e is the energy input due to 
the forcing. We thus have a new degree of universality for 
the 3/2 Kolmogorov law in two-dimensional turbulence, 
with respect to the class of equations (|l|) parameterized 
by A. From a dimensional point of view, (|^) implies a 
scaling exponent h — 1/3 for velocity increments and 
thus a scaling exponent 4/3 for ip, as in Navier-Stokes 
turbulence. From (^) one obtains the different predic- 
tions for the spectrum discussed above. In particular, 
the kinetic spectrum has the form 



(7) 



for any value of A, but with a Kolmogorov constant C\ 
which, in principle, can depend on A. A simple physical 
argument for this dependency is as follows. The scal- 
ing of the eddy turnover time depends on the scaling 
exponent. In the kinetic limit k 3> A one has the stan- 
dard Kolmogorov scaling r(fc) ~ e _1 / 3 fc~ 2 / 3 [n3[. On 
the other hand, in the potential limit k <C A, (Uj) gives 
r(fc) ~ (A/fc)^- 1 / 3 *;- 2 / 3 . Thus, for k < A, the efficence 
of energy transfer is reduced and one expects a larger 
value of the constant in (Q) . 

We have numerically investigated the inverse cascade 
in the potential energy regime by direct numerical simu- 
lations. In order to avoid complications induced by the 
crossover from the kinetic domain to the potential do- 
main, we study the system in the limit A — > oo. This is 
to be seen only as a formal procedure, equivalent to con- 
sidering wavenumbers much smaller than A, which physi- 
cally might be the case for magnetized plasma in presence 
of a strong magnetic field. Indeed this limit provides us 
with a model suitable for any A > fc/: because the en- 
ergy is transferred to large scale, the dominance of the 
potential term will be assured in all the inertial range. 
In the limit A ^> kt, rescaling the time t — > t/X 2 , one 
obtains the so-called asymptotic model M 



^ = J{V 2 ^^) + D + F 
for which the conserved quantities become 



(8) 



(9) 
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spectrum E\(k) at times 
(x) and t — 10"^ (+) and kinetic energy spec- 
trum Ek(k) at time t = 1CP 2 (*) averaged over 14 realiza- 
tions. The continuous line represents the dimensional pre- 
diction fc -11,/3 , the dashed line is the Kolmogorov spectrum 
E k (k) = C\e 2/3 k~ s/3 with constant C\ = 11. In the inset we 



plot the average energy flux at time t — 10 
e = 1.24. 



with the line 



We have integrated (||) with a standard pseudo- 
spectral code in a double periodic domain of size L = 2tt 
at resolution N — 512. The forcing is white in time in 
a narrow band of wavenumbers around kf — 160. The 
dissipative term in (j^) has the role of removing potential 
enstrophy at small scales and, as customary, it is numer- 
ically substituted by a hyperviscous term (of order 8 in 
our simulations). Time evolution is obtained by a stan- 
dard second-order Runge-Kutta scheme starting from a 
zero initial condition. The run is stopped at a given time 
T at which the energy containing scales are still much 
smaller than the computational box in order to avoid 
condensation effects H (see Figure [l]). All the results 
discussed in the following are taken after averaging over 
14 independent realizations. 

The limitation in the resolution (N = 512) is due to 
the discussed scaling of the characteristic time. Even 
with this moderate resolution, the ratio of the large scale 
characteristic time with the forcing scale time is about 
2000 and thus time evolution is very expensive (10 6 time 
steps for each realization) . In the case of Navier-Stokes 
turbulence (A = 0) this would correspond to an integra- 
tion covering about 5 decades of inertial range. 
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In Figure [l] we plot the potential energy spectrum 
E\(k) at two different times. The scaling exponent 
— 11/3 is clearly visible even if some accumulation at 
the largest mode is evident. This accumulation is not 
due to condensation as it is still well below the largest 
mode and it moves in time. We think that the existence 
of this "bump" is a genuine effect, probably due to the 
rapid growth of characteristic times and to the presence 
of intense vortices, as discussed below. The energy flux 
£ ~ 1.24 is estimated by the plateau of the energy flux 
shown in the inset. 
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FIG. 2. Third-order longitudinal structure function 
((5ml(Y)) 3 ) compensated with the dimensional prediction er 
at times t = 4x 10~ 3 (x) and t = 10~ 2 (+) averaged over 14 
independent realizations. The line represents the Kolmogorov 
law @. 



The "3/2" law for the longitudinal velocity structure 
function is shown in Figure g, also plotted at two differ- 
ent times. The compensation with the theoretical pre- 
diction (||) is remarkable, taking into account the limited 
resolution of our runs. As expected, the extension of the 
inertial range increases with time without changing small 
scale statistics. The oscillations observed at small scales 
are due to the contamination of the forcing. A similar 
effect was observed also in NS simulations. 

As discussed above, the fact that the "3/2" law is in- 
dependent on A (and thus the velocity scaling exponent 
has always the Kolmogorov value h = 1/3) do not im- 
ply that the statistics, and in particular the form of the 
pdf of velocity differences, is the same as for Navier- 
Stokes equation. For example, in Figure ^ we also plot 
the kinetic energy spectrum Ek(k) at final time t — 10~ 2 . 
The scaling exponent is compatible with the Kolmogorov 
value 5/3 as predicted by (M), but the Kolmogorov con- 
stant Ca — 11 is about two times that of Navier-Stokes 
Pj . A larger constant means a suppression of the energy 
flux which is a direct consequence of the dilatation of the 
characteristic times. 




FIG. 3. Grayscale plot of the stream function f/j at time 
t = 10 -2 . The characteristic scale of vortices £e ~ L/10 
corresponds to the peak of the energy spectrum in Figure [[]. 

A more significant difference with respect to Navier- 
Stokes inverse cascade is the presence of strong vortices, 
as shown in Figure 0. Vortices in CHM turbulence are 
injected at the forcing scale and they organize themselves 
to form a random pattern on the characteristic scale £e- 
This dynamical scale is associated to the peak of the 
spectrum of Figure H. Vortex dynamics slows down as 

— 8/3 

Ie increases as t(^) ~ . Thus, in the limit of large 
Reynolds number the system will end in a disordered pat- 
tern of quenched vortices forming a kind of "turbulent 
glass" . 

An important consequence of the presence of strong 
vortices is that the statistics of the velocity field strongly 
deviates from Gaussianity. In Figure || we plot the pdf of 
longitudinal and transverse velocity differences at three 
different scales within the inertial range. The effect of 
vortices is evident by the presence of large "wings" in 
the tails, in particular on the transverse velocity differ- 
ences which are more sensible to a rotating structure. 

In conclusion, we have shown that Kolmogorov "3/2" 
law for two-dimensional energy cascade in Charney- 
Hasegawa-Mima turbulence is independent on the value 
of the intrinsic scale A. Velocity statistics satisfies Kol- 
mogorov scaling with non- universal coefficients. In the 
asymptotic limit A — > oo the Kolmogorov constant is 
found to be about 2 times the Navier-Stokes case [ p)| ]. 
Strong coherent vortices are found to emerge at the forc- 
ing scale and aggregate to form a pattern of quenched 
vortices at large scale. As a consequence of the pres- 
ence of vortices, strong deviations from Gaussianity are 
observed in the velocity field. 
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FIG. 4. Probability density functions of longitudinal (a) 
and transverse (b) velocity differences at separations I — 0.05 
(+), I = 0.1 (x) and i = 0.2 (*) at time t = 10~ 2 . Dashed 
line represents Gaussian distribution. 
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